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We investigate the effect of imperfections in realistic detectors upon the problem of quantum state 
and parameter estimation by continuous monitoring of an open quantum system. Specifically, we 
have reexamined the system of a two-level atom with an unknown Rabi frequency introduced by 
Gambetta and Wiseman [Phys. Rev. A 64, 042105 (2001)]. We consider only direct photodetection 
and use the realistic quantum trajectory theory reported by Warszawski, Wiseman, and Mabuchi 
[Phys. Rev. A 65, 023802 (2002)]. The most significant effect comes from a finite bandwidth, 
corresponding to an uncertainty in the response time of the photodiode. Unless the bandwidth is 
significantly greater than the Rabi frequency, the observer's ability to obtain information about the 
unknown Rabi frequency, and about the state of the atom, is severely compromised. This result 
has implications for quantum control in the presence of unknown parameters for realistic detectors, 
and even for ideal detectors, as it implies that most of the information in the measurement record 
is contained in the precise timing of the detections. 
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I. INTRODUCTION 

Quantum parameter estimation originated with the 
work by Helevo [1] and Helstrom [2] and is usually for- 
mulated as follows. A known quantum state undergoes 
some evolution. This evolution is usually unitary but 
need not be [3], and is parameterized by one or more 
unknown parameters. The goal is to estimate these pa- 
rameters by making a measurement on the final state. In 
a recent paper, two of us [4] , following Mabuchi [5] , con- 
sidered the situation in which an open quantum system 
parameterized by an unknown dynamical parameter is 
estimated through continuous-in-time measurements of 
the environment. Ref. [4] also focussed on the condi- 
tioned state and how the lack of information about the 
parameter affects its purity. These two papers are part 
of a small but growing body of work including [6] and 
[7] (and references therein) that address the problem of 
reducing classical uncertainties in the evolution of quan- 
tum systems by performing continuous-in-time measure- 
ments. 

In Ref. [4], the system is taken to be a damped, classi- 
cally driven two-level atom (TLA) with a driving Hamil- 
tonian (oc a x in a rotating frame) with an uncertain driv- 
ing strength O. The physical scenario the authors give 
for this is as follows. An atom is placed in a classical 
standing wave, thus experiencing a Rabi frequency of 

& = fi max sin(fcr), (1.1) 

where k is the wavevector for the classical field and x 
is the position of centre of mass of the atom. Gam- 
betta and Wiseman then assume that the TLA is equally 
likely to be placed at any position in the standing wave. 
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Hence, a range of — Sl m ax < ^ < fimax exists for the 
possible driving strength. It is also assumed that the 
TLA has zero translational motion. This can be justified 
by assuming that the TLA is fixed by a radio frequency 
Paul trap or similar confining mechanism [5]. Although 
it is shown that continual efficient measurement upon 
the output electromagnetic field of a TLA leads to in- 
formation gain about the driving strength, the amount 
of information gained depends strongly on the measure- 
ment scheme used. The purity of the conditional state 
also depends strongly on the measurement scheme, but 
the performance of schemes by these different criteria are 
not well correlated at all [4] . 

In this paper, we consider only one sort of detection: 
direct detection, but we consider imperfect detection. 
This allows us to address the following two questions. 
First, what (qualitatively) are the aspects of the ideal 
measurement record that contain the information about 
the system state and the unknown parameter? Second, 
in practice, how badly would the amount of information 
gained be affected by detectors with realistic imperfec- 
tions. To treat realistic detection we use the theory de- 
veloped in Rcfs. [8, 9, 10]. This generalized the theory of 
quantum trajectories [11, 12, 13, 14, 15, 16] by determin- 
ing the state of the system conditioned upon the output 
of a detector that has dead-time (rdd)j dark counts at 
rate ydk, inefficiency rj, and finite bandwidth 7 r . A dia- 
gram of this detector is displayed in Fig. 1. 

The reciprocal of the bandwidth 7" 1 sets the scale for 
the uncertainty in the delay between the absorption of 
photon (which is synonymous with its detection in an 
ideal detector) and the avalanche of charge which reg- 
isters this at a macroscopic level in a realistic detector. 
We find that a finite bandwidth has the greatest effect 
on our parameter estimation problem. This tells us that 
it is the timing between detections, not the steady state 
detection rate, that contains the majority of the infor- 
mation gained. Unless the bandwidth of the detector is 
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REALISTIC PHOTON COUNTER 




FIG. 1: A schematic illustrating realistic direct detection. 
The atom is placed at the focus of a parabolic mirror such 
that all the fluorescence enters the realistic detector. Realistic 
photon counting is performed by an avalanche photodiode. 
The quantum efficiency r] of the diode is represented by the 
the beam-splitter (BS). Single arrowheads within the realistic 
photodetector indicate Poisson processes. The microscopic 
detector states are indicated by a 0, 1, and 2. 



substantially greater than f2 max , the rate of gain of infor- 
mation about f2 is much reduced, and consequently the 
purity of the system is also. Like the work of Ref. [4] this 
has important implications for the robustness of quantum 
feedback control (see for example Refs. [17, 18]). 

This paper is structured as follows. In Sec. II we be- 
gin with some numerical results to motivate the theo- 
retical developments in the remainder of the paper. In 
Sec. Ill we briefly review the theory of parameter esti- 
mation based on linear quantum trajectories in Ref. [4], 
and show how it can be applied to realistic detection. In 
Sec. IV we use this theory to investigate the evolution of 
the conditional state for an initially unknown parameter, 
and in Sec. V we investigate the rate of information gain 
about the parameter. We conclude with a discussion in 
Sec. VI. 



II. CONDITIONAL PROBABILITY 
DISTRIBUTION FOR DRIVING STRENGTH 

Before providing a theoretical basis, some numerical 
results are provided to give the reader an insight into 
the effects of realistic detection on parameter estimation. 
With an example in mind, the theory contained in the 
following section should also prove to be more transpar- 
ent. 

We begin by looking at how an observer's probabil- 
ity distribution for the driving strength evolves given a 
single measurement record (trajectory). Details of the 
numerical techniques employed to do this will be given 
at the end of the next section. The distribution is de- 
noted by P(f2|I[o j t)), where I[o,t) represents the measure- 
ment record. Obviously, I[ 0i t) refers to different events for 
ideal and realistic detection: it is a record of the times 
of photodetections and of avalanches, respectively. The 
initial distribution can be found from Eq. (1.1) with x 
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FIG. 2: Plots (A) and (B) show how the probability distribu- 
tion for Q typically evolves while continuous direct detection 
is performed on the output of the quantum system (a TLA 
inside an optical cavity). Plot (A) is for ideal detection, while 
(B) is for realistic detection. The parameters used for the pho- 
ton counter were r\ = 100%, and 7 r = 7F, r d d = 0, 7dk = 0. 
The value of Q used to generate I[o,t) was fitruc = 4r. The 
maximum allowed value for Q is f2 max = 10r. 



taken to have a flat distribution. The result is 
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(2.1) 



The unconditioned evolution of the atom (which is inde- 
pendent of the scheme, or quality, of detection) is given 
by the master equation (ME) 



P(*) = -y[o*,p(<)] 



TV[a]p(t) = C a p{t). (2.2) 



Here fl is assumed known, T is the spontaneous emission 
rate (which we will set equal to unity), a is the atomic 
lowering operator, and a x = a + a' . The superoperator 
T> gives the damping of the system into the environment 
and is defined as [19] 



T)[a]p = apa) — h {a^ap + pa) a] 



(2.3) 



A comparison of ideal and realistic detection is given in 
Fig. 2. In this figure Qtrue = 4r, but the observer does 
not know this, but rather has knowledge described by the 
prior distribution (2.1) with O max = 10T. As time goes 
on, the detection record reveals more about f2true, and 
this probability is updated in a Baycsian manner. Note 
that for both ideal and realistic detectors, the distribu- 
tion always remains symmetric. This is because the sign 
of Q cannot be determined from direct detection. Pho- 
toemissions depend only upon the mean of a z , and that is 
independent of whether the TLA state rotates clockwise 
or counterclockwise around the a x axis. 

In Fig. 2, we have chosen a realistic detector whose 
only imperfection is a finite bandwidth equal to 7r (see 
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figure caption). It can be seen that the degradation of 
the quality of measurement is quite large even though our 
'realistic' detector is well beyond current technology. If 
we choose the realistic parameters for an APD (avalanche 
photodioide) used in Refs. [8, 10] we obtain results that 
are so far from those of ideal detection that a comparison 
lacks interest. 

The fact that the chosen detector parameters give re- 
sults that are still poor compared to ideal measurement 
is a little surprising. Very few photodetections are missed 
under realistic detection with a finite j r being the only 
imperfection. It illustrates that the rate of information 
gain is very sensitive to an uncertainty in the state of the 
TLA after a jump. 

The plot for realistic detection, Fig. 2 (B), also shows 
that it is difficult to distinguish between Sltrue and an in- 
teger multiple of fitruej say mf2 t ruc (for this figure this oc- 
curs only for |f2t rU o| and 2|f2 t mc|)- By fitrue we mean the 
actual driving strength that exists, but which is initially 
unknown to the observer. This difficulty exists because 
the evolution for these two proposed values of the driving 
strength will place the TLA in the excited state at the 
same time every integer oscillation (at the smaller fre- 
quency). Hence, for certain times there will be a strong 
correlation between the probabilities of decay. Obviously, 
a photoemission can still occur when a driving strength 
of 2|f2truc| would place the TLA very near ground state, 
thus ruling out that possible value. Whether these de- 
cays that are discriminatory in this sense are detected or 
not depends on the choice (and duration) of the quantum 
trajectory in question. In Fig. 2 (B) a small but, persis- 
tent (at least up to t = 50) probability for = 2|f2 truo | 
exists. 

It is worth mentioning that the plots in Fig. 2 are based 
on the 'same' photo-absorptions in the ideal and realistic 
cases. This is possible using the simulation technique 
of Refs. [10], in which the realistic quantum trajectory 
can be assured to be consistent with the ideal quantum 
trajectory. 



III. PARAMETER AND STATE ESTIMATION 
THEORY 



We now briefly present the theory of parameter and 
state estimation in the context that is relevant to this 
paper. The specific problem that we wish to address 
involves calculating the new probability distribution for 

and also obtaining the best estimate of the TLA state, 
based on the measurement record. Bayesian statistics are 
required, with the new distribution, based on L t ), given 
by [20] 



The best estimate of the TLA state is found from [4] 



p(i m \n)p (n) 
J p(i [0 , t) \n)p (n)dn- 



pi(t) = 



J pi,n(t)Pn(n)dn 

J p(i M \n)p (n)dn : 



(3.2) 



where pi t n{t) is the unnormalised state of the TLA given 
that there was a measurement record I\o,t) and that the 
driving strength was fi. 

In this paper we leave the state unnormalised even af- 
ter jumps, so that 



(3.3) 



meaning that P(I[ 0t ) \Q) is actually equal to the norm of 
the TLA state, when it is evolved with Q under I[o,t)- 
Ideally I[o,t) would be an experimental record, but for 
convenience we simulate it using fitmc and then 'forget- 
ting' it. The integrals in Eqs. (3.1)-(3.2) require the sim- 
ulation of the evolution for the entire range of f2. This 
means that the process is very time consuming as many 
trajectories (for the different fi's) have to be run in order 
to obtain a single trajectory for P(f2|I[o it )) or pi(t). 

In order to treat realistic detection, it is necessary to 
describe not just the state of the TLA, but also the state 
of a detector (treated classically) which is correlated with 
the atomic state [8, 9]. We do this by introducing a 
set of three operators pi, s (t) for this supcrsystem, where 
here s = 0,1,2 represents the three classical detector 
states. These correspond to the ready state (s = 0), 
the avalanching state (s = 1) and the resetting state 
(s = 2). Refer to Fig. 1 for details. The state of the TLA 
is defined as pi(t) = pi,o{t) + pi,i(t) + pi^{t)- Equation 
(3.2) is trivially generalized to give the best estimate of 
the supersystem state pi tS (t): 



pi At) = 



J Pi,a,s(t)Po(tydn 
f Pp m \(l)P (Sl)dSl' 



(3.4) 



(3.1) 



The three equations, Eqs. (3.1)-(3.2) and Eq. (3.4), 
are unusable in their present form as P(I[ 0jt )|f2) and the 
norm of pi,n(t) (which is equal to P(L 0ii ) |fi)) will be- 
come extremely small as t increases. This will be true 
even for Q = Otrue; as can be seen by considering the 
probability of a string of measurement results occurring. 
The probability of a jump in a given infinitesimal inter- 
val is never greater than Tdt, so that the norm of the 
state after m jumps would be less than (Tdt) m (since the 
no-jump evolution also decreases the norm of the state). 
In simulations, a time step of dt = 10 _4 r _1 is typical. 
Thus, numerical error introduced by computer round-off 
necessitates a different approach which we now explain. 



A. Linear Quantum Trajectories 

In linear quantum trajectories [16], the state is multi- 
plied by a predetermined factor after every measurement 
event. The normalisation constants arc chosen so that 
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the norm of the states corresponding to the most likely 
values of ft stay relatively close to unity. Each possible 
measurement result, r, has an 'ostensible' [16] probabil- 
ity, A(r) associated with it, that we use as the normalis- 
ing factor. Therefore, we have 



Pi,n,s{t) 



Pi,n, s (t) 
A(I[o,t))' 



(3.5) 



where pi t u,s is the supersystem state when linear quan- 
tum trajectories are used and A(Iro,t)) is the ostensible 
probability for getting I[o,t)) 



A(I [0 , t) )=A(r fe )A(r fc _ 1 )...A(r 1 ). 



(3.6) 



Here, we have assumed that the time interval [0, t) has 
been divided into a very large number, k = t/dt, of dis- 
crete measurement times. The actual probability of get- 
ting I[o,t); given ft is 



(3.7) 



where the trace is over the TLA and the summation is 
over the detector states. 

This allows Eq. (3.1) and Eq. (3.4) to be written as [4] 



and 



PW[o,t)) 



Pl,s(t) 



E s Tr[pi,o,,(t)]Po(ft) 
/E.'fr[frA.(*)]fl>(fi)dfi 



J ft,n, a ft)fo(n)dn 

JEs^[pi,nAt)}P a (n)dn- 



(3.8) 



(3.9) 



Since Trfp^n,^)] i s °f order unity for the most likely 
ft's the problems associated with using Tr[pi^ s (t)] have 
been avoided. The best estimate of the system state 
allows the conditional evolution of a single stochastic 
record to be investigated, while P(ft|I[ 0jt )) will allow cal- 
culation of the information gain (see following). 

In this paper, we choose the ostensible probabilities to 

be 



A(l) = edt, 



(3.10) 



where r = 1 indicates the detection of a photon (ideal 
detector) or an avalanche (realistic detector). No detec- 
tion (r = 0) has an ostensible probability A(0) = 1 — edt. 
With these ostensible probabilities, the linear quantum 
trajectory equations for the supersystem describing real- 
istic detection [8, 9] become 

dpo(t) = dt(£n-7dk-f/r i 7[5-] + e)po(t) 

- dN{t)p {t)+dN{t){t - r dd )p 2 {t) (3.11) 

dptit) = dt(C n - ll + e)p 1 {t)+dt{i 1 rj[a}+j dk ) 

x poi^-dJVit)^ (3.12) 

dp 2 (t) = dt(jC Q + e)p 2 {t)+dAf(t)p 1 (t)/e 

-dN{t){t-T dd )p 2 {t). (3.13) 



Here dN{t) = 0,1, the increment in the number of 
avalanches, is equivalent to the measurement result r in 
the time interval [t, t + dt). We maintain the symbol I[ 0j t) 
for the record of avalanches in the interval [0, i), but for 
clarity we have omitted the subscripts I and ft for the 
p(t)'s. The superoperator J is defined by [19] 



J[a]p = apa) . 
For simulation purposes we take 



n 



2^-ue + r 2 ' 



(3.14) 



(3.15) 



the expected detection rate in steady state (ignoring any 
dead time) when ft is known. This is allowed as the 
results of the simulation are independent of e, thus we 
are free to choose any value. Note that renormalisation 
of p(t) is occurring even when the detector is resetting 
and the probability of an avalanche is strictly zero. This 
illustrates the flexibility of the ostensible probabilities, 
where the only concern is maintaining a norm of order 
unity. 



IV. CONDITIONAL STATE EVOLUTION 

Now that we have explained how the numerical results 
of Sec. II were obtained, we return to these simulations 
to investigate the conditional state evolution in a typi- 
cal realistic quantum trajectory. In Fig. 3 we see how 
z(t) (the mean of a z ) behaves under a variety of cir- 
cumstances. Specifically, we consider known |fttrue| and 
ideal detection (thin dashed line), unknown |fttrue| (that 
is, an initial distribution given by Eq. (2.1)) and ideal 
detection (thin solid line), known | f2truc | and realistic de- 
tection (thick dashed line) and, finally, unknown fttrucl 
and realistic detection (thick solid line). The reason why 
a comparison is made between the trajectories for a TLA 
randomly placed in the standing wave of the field and 
trajectories of known |f2true|> but unknown sign, is that 
at t = oo the sign of fttruc will still not be known, as 
discussed in Sec. II. This allows us to reveal how far 
the initially unknown |fttrue| trajectories are from their 
optimal form, given the constraints of direct detection. 

The conditioned state corresponding to a known |fttrue| 
is formed from 



Piit) = (pi,n true (t)+Pi,-n true (t))/2. 



(4.1) 



That is, it is an equal mixture of the two conditioned 
states obtained with a driving strength of +f2true an d 
— fttruc- Given that we take the TLA as starting in the 
ground state, these two conditioned states will have the 
same z(t), but opposite sign for y(t) resulting in y(t) = 
for the mixture as initially Pq(Q) is an even mix of posi- 
tive and negative ft. If the initial state of the TLA (which 
is presumed to be known) gives a finite value for y(t = 0) 
then the time of the first jump will give a very small 
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FIG. 3: Two segments of evolution of z(t) in a typical trajec- 
tory are shown. The legend is as follows: thin-lines are for 
ideal detection, thick-lines are for realistic detection, dashed 
is for known |Otrue| and solid is for unknown |r2 t ruo|- Thus, to 
be explicit, the thin dashed line is for |S!truc| known and ideal 
detection, the thin solid line is for |f2true| unknown and ideal 
detection, the thick dashed line is for |f2true| known and re- 
alistic detection and, finally, the thick solid line is for |fi trU e| 
unknown and realistic detection. The detector and system 
parameters (and time units) are as for Fig. 2. 



amount of information about f2 trU c- This is because an 
early jump would tend to support the conclusion that the 
TLA had been rotated on the Bloch sphere in the direc- 
tion that takes it towards the excited state first, while a 
later jump suggests that the TLA was rotated in the op- 
posite direction. However, the initial state is taken to be 
the ground state in our simulations. Since for direction 
detection with no detuning x(t) =0 [11, 19], we have 



pi(t) = (I + z{t)a z )/2, 



(4.2) 



where z{t) could be obtained from either the +f2t ruc or 
— f^true trajectory. It should now be clear why only z(t) 
is plotted. In fact, the purity can also be inferred from 
the z(t) plot as pi(t) = (1 + z 2 (t))/2. 

As for Fig. 2, the realistic trajectory is forced to be 
consistent with the ideal trajectory. This is evidenced 
by the jump times being correlated. Because of the de- 
tail of these plots we show two portions of the evolution 
only, illustrating how improved knowledge over time of 
the dynamical parameter, fl, allows closer adhesion to 
the maximum knowledge trajectory (thin dashed line for 
ideal detection and thick dashed line for realistic detec- 
tion). 

In the first 10F -1 units of time the frequency of oscil- 
lation of z(t) for the realistic detection state (solid black 
line) is poorly defined. One can also see that the fre- 
quency is faster than the true frequency, due to Po(^) 
being peaked at ±f2 max [4]. The small amplitude of os- 
cillation is also due to the lack of knowledge about fl. 
The range of possible frequencies ensures that when one 



value of would place the TLA in the excited state, an- 
other would place it near the ground state, thus averaging 
the amplitude of oscillation of z(t) to below unity. 

For ideal detection (thin solid line), the frequency of 
oscillation is also shifting about, but the amplitude of 
the oscillations have already noticeably increased in size 
by t = 10r _1 . Another feature is that jumps take z(t) 
to —1 (a pure state) for ideal detection even if r^trucl 
is not known. This is because there is no time delay in 
conveying the decay of the TLA to the observer. This is 
not true for realistic detection where the jump is slightly 
delayed and the TLA will have rotated out of the ground 
state by some unknown amount. 

After 40r~ 1 time units, the ideal detection observer 
has gained enough information about |f2true| that z(t) 
matches well with the 'ideal' z(t) (Fig. 3 (B)). The re- 
alistic observer's z(t) is much closer to the ideal than 
at the start of the trajectory, but some uncertainty still 
obviously exists. 



V. INFORMATION GAIN 

We now move away from conditional dynamics and 
look at how the ensemble averaged gain of information 
concerning fl is affected by realistic detection. As in [4], 
the measure used to quantify the quality of information 
gained about fl is AIj. It is given, in bits, by 

Ah = f dfiP(fi|I [0 , t) )log 2 P(fi|I [0>t )) 

dnp {n)\og 2 p (n). (5.1) 

As the distribution for f2 becomes more sharply peaked 
around |f2 t rue|j All will increase. Due to the stochas- 
tic nature of I[o.t), an ensemble of trajectories are run in 
order to obtain the average, AI, of Ah that is charac- 
teristic of the measurement scheme. Each member of the 
ensemble is formed by picking an O trU e randomly (but ac- 
cording to the distribution Eq. (2.1)) and simulating the 
evolution based on a stochastic trajectory corresponding 
to that f2truc- A single trajectory for each member of an 
ensemble of fame's is sufficient as the averaging of the 
stochasticity is done over the multitude of fame's- 

A major difference between parameter estimation for 
ideal direct detection and realistic direct detection is 
found when f^truc 3> L. We investigate this regime by 
varying the value of f2 max used in the simulations. This 
works as the most likely values for f2truc come from close 
to f^max- By increasing f2 max we are effectively increas- 
ing the values of £ltmc used. This method corresponds 
physically to the TLA being placed randomly in standing 
waves of different amplitude. 

In Fig. 4 (A) the information gain, AJ, averaged over 
an ensemble size of 1000 (that is, 1000 different roue's), 
is shown for ideal detection. The three plots are for 
^max = 5r, 10r, 20r. In Fig. 4 (B) the results are shown 
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FIG. 4: In plot (A) the ensemble averaged information gain 
for ideal detection is shown for three different values of f2 max . 
In plot (B) the ensemble averaged (over 100 trajectories) in- 
formation gain for ideal detection is shown for the same values 
of fimax- In plot (C) the difference between ideal and realistic 
detection is shown. 



for the same values of r2 max when the monitoring is done 
with realistic detection. Error bars are included as the 
ensemble size in this case is only 100, due to the greater 
computational intensity. As f2 max doubles, the informa- 
tion gain curves for ideal detection are separated by ap- 
proximately 1 bit at large t. This does not indicate that 
the detection scheme is more effective at large f2 max , 
rather, it is an artifact of the initial information con- 
tent of the distribution Po(f2). This can easily be shown 
analytically. Consider the integral 



dfl- 



1 



log 2 



1 



2 

max 



n 2 



(5.2) 



If the change Q' msix — > mf2 max is made, a change of in- 
tegration variable to CI' = mfl allows the information of 
the new distribution to be expressed in terms of the old 
one. The result is 



Ah 



I,(mO„ 



AL 



i,(n„ 



log 2 TO. 



(5.3) 



Thus, every time ri max is doubled, the initial informa- 
tion decreases by 1 bit. Since the initial information is 
being subtracted from the information at time t, we can 
conclude that the width of the P(il) distributions, for 
different f2 m ax, in the long time limit are roughly equal 
for ideal detection. 

This is qualitatively different from the case of realis- 
tic detection where the information gain for f2 max = 5I 1 
and 10r are roughly the same and f2 m ax = 20r sees a 
significant decrease. Taking into account the fact that as 
!] max increases less information is being subtracted away 
from the information at time t, it is clear that realistic 
detection is becoming considerably worse at distinguish- 
ing between candidates for Ot rue . This is more clearly 



illustrated in Fig. 4 (C) which shows the difference be- 
tween ideal and realistic detection information gain for 
the three values of fi max . We now try to understand 
these trends by discussing what features of a measure- 
ment record allow the information to be accrued. 

The most simple (and naive) way to determine the ac- 
tual driving strength is to analyse the mean steady-state 
photon counting rate (I) ss . For ideal direct detection this 
is 



TQ 2 



20? + r 2 ' 



(5.4) 



By using Eq. (5.4) the mean flux can be matched with a 
particular £1, given that T is known. However, the flux 
asymptotes to r/2 at large fi, essentially making one 
value of Q indistinguishable from another. Because ideal 
detection still gives information at large f2 max we deduce 
that the specific times at which detections are occurring 
also play a role in determining f^true- 

To consider this effect in more detail we look at the 
waiting time distribution for photon detections. In per- 
fect direct detection [21] 



w(t) 



o 2 



n 2 - (r/2) 



^-^sin^ip 2 -^) 2 ] 1 /^}, 

(5.5) 

where w(t) is the probability density of a photon be- 
ing emitted by the TLA at a time r after the previous 
emission. The waiting time, w(t), is zero at r = as 
the TLA in the ground state after a detection. A plot 
of this function for Q = ST. with T = 1 is given in 
Fig. 5 (A). The point that is meant to be taken from 
this graph is that there are oscillations that make certain 
emission times more likely (when the TLA is in the ex- 
cited state) and other times prohibited (when it is in the 
ground state), given a certain driving strength. These 
oscillations, which take w(t) to zero, are present for all 
51 > r/2, which includes the regime of interest here. This 
allows the narrowing of the distribution for O, based on 
a measurement record, despite the fact that very little 
information can be obtained from the mean flux rate. 

The waiting time distribution for realistic detection is 
much more complex. In this case, r would be the waiting 
time between avalanches of the photodiodc, rather than 
emissions of the TLA. One of the difficulties lies in the 
fact that there is not a definite state of the atom after an 
avalanche, due to the finite response time. Thus, initial 
conditions for the TLA must be chosen as being charac- 
teristic of the post-avalanche state. We chose them to be 
the solution to the ME (2.2) averaged over the time for 
the avalanche to mature. That is. 



Paval. = 7r 



-7r* 



P (t)dt, 



(5.6) 



where p{t) is the solution to the ME given that at t = 
the TLA is in the ground state. Of course, this will only 
be true if the dark counts arc ignored, which we will do. 
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The integral in Eq. (5.6) can be carried out analytically 
with x — and the y and z components of the state being 



yaval. 



^aval. — 



20(1 + 7r 



(l + 3 7r + 2 7 2 + 2S7 2 )' 
(2 7r + l)(l + 7r ) 
(1 + 3 7r + 2 7r 2 + 2S1 2 ) ' 



(5.7) 
(5.8) 



where we have set T = 1. 

If we take the simplest case where Tdd = and ^dk = 0, 
then we can use these as the initial conditions for the 
ready state of the detector. (The case where tm ^ is 
only a little more complicated, because the solutions in 
Eqs. (5.7)-(5.8) can be evolved forward a time r^d with 
the ME to give the state of the TLA when the detector 
becomes ready to detect photons once more.) A series of 
6 linear coupled differential equations must be solved in 
order to determine 7r Tr[pi (r)] = w(r) ava i. , which is equal 
to the probability of an avalanche at time t (the waiting 
time). Here, Pi(t) is the solution to the unnormaliscd 
evolution equations given in Eqs. (3.11) and (3.12) with 
e = and dAf = (no avalanches) from Tdd to r. 

The 6 real variables that are necessary can be expressed 
as a vector r= (yo, zo, Po, yi, Z\, andPi) T , where the the 
subscripts refer to the detector state and the P's are the 
probability of occupation. The equations of motion for 
this vector is dr/dt = At where A is the following matrix 
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-1 



~7r 



1 



(5.9) 
The 



whereas before we are assuming 7 dk = and T 
variable of interest is then io(r) ava i = 7r Pi(r) 

Plots of iy(r)avai. with SI — 5r, 10r, and 20r are given 
in Fig. 5 (B), (C) and (D). Most notable is that although 
there are oscillations, they do not take iiWaiXi") to zero. 
The random response time of the photodetector is wash- 
ing out the peaks and troughs of the ideal detection wait- 
ing time distribution, making it much more difficult to 
distinguish between Si's on the basis of the specific times 
of the avalanches. This effect becomes more pronounced 
as S7 increases. The relevant characteristic time scale is 
the response time of the detector. For SI > 7r , the TLA 
may have rotated out of the ground state significantly 
by the time the avalanche occurs. If the state of the 
TLA is smeared out and information contained in the 
higher order moments of the waiting distribution is lost, 
then there is no way to determine S7 tm e| when it is large. 
This is the reason for the drop-off of the information gain 
for realistic detection that is seen in Fig. 4. 

Although we have done our simulations with a dead 
time Tdd = 0, the finite response bandwidth leads to an 
effective dead time of r )~ 1 , which leads to an effective 
drop in the efficiency of the detector. To rule inefficiency 
as the cause of the loss of structure in w(t) as S7 increases 




f 04 





5 10 

Time.x [r- 1 ] 



FIG. 5: In plot (A) the waiting time distribution for ideal 
detection is shown for SI = 5I\ In plots (B), (C) and (D) 
the waiting time between avalanches for realistic detection is 
shown for the values of Q = 5r, S7 = 10r and Q = 20I\ The 
photon counter parameters are as for Fig. 2. 





5 10 
Time,-! [r- 1 ] 



FIG. 6: In plot (A) the waiting time is plotted for Q = 5r 
and 7] = 100%. In plots (B), (C), and (D) the waiting time is 
plotted for S7 = 5r , S7 = lor and fi = 20r and rj = 70%. 



we consider the waiting time distribution for a detector 
of efficiency rj. Using the formula given in Ref. [21] we 
have calculated w(t) with ft = 5T , 10T. and 20r. Plots 
are given in Fig. 6 (B), (C) and (D) (with 77 = 70%) while 
Fig. 6 (A) is a reproduction of Fig. 5 (A) (S7 = 5r and 
rj = 100%). Here it is observed that the smearing out of 
w(t) (for detectors with only an inefficiency) is due only 
to this inefficiency. That is, it is independent of S7. Thus 
the smearing out the waiting time distribution as a result 
of increasing SI for a realistic detector is brought about 
by the finite bandwidth, 7r . 
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VI. DISCUSSION 

We have investigated the effect of imperfections in re- 
alistic detectors upon the problem of quantum state and 
parameter estimation by continuous monitoring of an 
open quantum system. Specifically, we have re-examined 
the system of a two-level atom (TLA) with an unknown 
driving strength introduced in Rcf. [4] . Considering only 
direct photodetection, we find that the most significant 
effect comes from a finite bandwidth of the detector. This 
corresponds to a randomness in the response time of the 
photodiode (the time from photon absorption to the time 
when the avalanche reaches a pre-set threshold). 

We find that unless the bandwidth is significantly 
greater than the Rabi frequency, the observer's ability to 
obtain information about the unknown Rabi frequency, 
and about the state of the atom, is severely compromised. 
This is because the waiting time distribution between 
avalanches is smeared out by the bandwidth, losing the 
Rabi oscillations that characterize it for an ideal detector 
(see Fig. 5). This result shows that the Bayesian update 
method in Ref. [4] implicitly made use of the structure 
of the waiting time distribution, rather than simply the 
mean detection rate. This implies that even for an ideal 
detector this parameter estimation problem is extremely 
nonlinear, and requires the full power of Bayesian prob- 
ability theory. Our results thus have implications for 
quantum control [17] in the presence of unknown param- 
eters even for ideal detectors. 

A final comment on the question of the physicality of 
the scenario that we have constructed needs to be made. 
In order to measure the emissions of an atom with an 
efficiency anywhere close to unity it would be necessary 
to place the atom in a micro cavity, with the cavity mode 



being heavily damped so that the atom can be regarded 
as emitting predominantly into the cavity output beam 
[22, 23]. However, in this paper we have considered a 
situation in which the position of the TLA is random, 
but fixed, across one wavelength of a standing wave field. 
This has the ramification that the effective decay rate is 
also position-dependent and hence also uncertain, scaling 
like the driving strength squared. Thus, in reality, there 
would be uncertainty in the two dynamical parameters 
f2 and T. In principle, with ideal detection, both these 
parameters could be found (apart from the sign of ft) 
given a long enough measurement record, as the waiting 
time distribution for ideal detection will be different for 
different positions of the TLA in the standing wave. It 
would be interesting to see how realistic direct detection 
affects the rate gain of information about £1 and F. 

Another question which would be interesting to answer 
is, given that homodyne detection of the y quadrature in 
[4] was shown to provide the greatest gain in information, 
does this result also apply when realistic detection affects 
are taken into account. Lastly, an interesting idea would 
be to investigate this scenario with a frequency-filtered 
measurement scheme, such as in Rcf. [24]. This would re- 
quire extending the measuring apparatus to include many 
cavities (one for each frequency) with each cavity having 
its own realistic detector. Either of these schemes would 
require far greater numerical resources than the (not in- 
significant) resources used here. 
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